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Landau damping in weakly interacting Bose gases is inves- 
tigated by means of perturbation theory. Our approach points 
out the crucial role played by Bose-Einstem condensation and 
yields an explicit expression for the decay rate of elementary 
excitations in both uniform and non uniform gases. System- 
atic results are derived for the phonon width in homogeneous 
gases interacting with repulsive forces. Special attention is 
given to the low and high temperature regimes. 
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The collective excitations of gases of alkali atoms con- 
fined in magnetic traps have been the object of recent 
experimental measurements |Q. The observed frequen- 
cies very well agree with the theoretical predictions of 
mean field theory confirming the validity of the Bo- 
goliubov scheme, later extended by Gross and Pitaevskii 
to the non homogeneous case, for describing the dynamic 
behaviour of dilute Bose gases. Conversely, the damping 
mechanisms associated with such excitations are not well 
understood and still represent a challenging question for 
theoretical investigation. 

The damping of collective modes can have various ori- 
gins. At T — O it can arise because of decay into two 
or more excitations with lower energy. This mechanism, 
which is well understood in uniform Bose superfluids || , 
is not active for the lowest modes of a trapped gas, be- 
cause of discretization of levels. Damping can be due 
to other non linear phenomena of classical or quantum 
nature Q. Experiments with magnetic traps however 
point out the occurrence of damping in an almost lin- 
ear regime where such effects should be negligible. At 
finite temperature damping can be due to collisional ef- 
fects of dissipative type. Thcsc effects are expected to 
be important at high temperature and density. Anothcr 
mechanism, holding in the collisionless regime, is Lan- 
dau damping which oecurs when the collective excita- 
tions can decay due to coupling with transitions associ- 
ated with other elementary excitations and occurring at 
the same frequency. Landau damping is not associated 
with thcrmalization processes and can be well deseribed 
in the framework of mean field theory (see, for exam- 
ple, U §§28-29, 34 and g). In a Bose gas it arises only 
at finite temperature. Its possible relevance to explain 
the experimental data in trapped Bose gases has been 
recently proposed by Liu and Schieve @ . 



The purpose of this lcttcr is to develop the microseopic 
formalism of Landau damping in a dilute Bose gas. The 
aim is twofold. On the one hand we provide a system- 
atic deseription of Landau damping in homogeneous Sys- 
tems covering the regimes of low temperature, first in- 
vestigated by Hohenberg and Martin ||, and the one 
at higher temperatures first investigated by Szepfalusy 
and Kondor H|. On the other hand we provide a for- 
malism suitable for the caleulation of Landau damping 
in trapped gases where the classification of elementary 
excitations significantly differs from the one of uniform 
systems. Our formalism is based on perturbation theory 
and points out the crucial role played by Bose-Einstein 
condensation. Previous approaches were based on kinetic 
theory for superfluids H or on the use of Green's funetion 
techniques B. 

Let E be the energy of the system associated with the 
occurrence of a classical oscillation of frequency oj in- 
duced, for example, by some external drive as happens 
in the experiments of ref. jl| . By classical oscillation we 
mean that the number of quanta of oscillation is very 
large. Due to interaetion effects, the thermal component 
of the gas can either absorb or emit quanta of this oscil- 
lation, thereby giving rise to the following expression for 
the energy loss 



E os = -Tiuj{W {a) ~ w {e) ) 



(1) 



In the above equation W^ a ' and are, respectively, 

the probabilities of absorption and emission of the cor- 
responding quantum huj given, in perturbation theory, 
by 



W = tt 



I (k | V„ 



(2) 



where the matrbc element is associated with a transition 
in which the i-th excitation, availablc in the system be- 
cause of thermal aetivation, is transformed into the fc-th 
one and Ek — Ei + huj in the case of absorption (a), and 
Ek = Ei~ hu in the case of emission (e) . 

In second quantization the interaetion term can be 
written in the following way 



V int = - I dr ^ ^ ipip , 



(3) 



where the coupling constant g is related to the s- 
wave scattering length through the usual formula g = 



1 



Anh 2 a/m. The field operator ip can be written as ifi = 
ipo + 5ip where ipo is the order parameter at equilibrium 
while 6ip characterizes its fluctuations which can be ex- 
pressed in terms of the annihilation (a) and creation (cr) 
operators relative to the elementary modes of the system 
in the usual form 



(r)at 



(4) 



and u, v can be determined, for example, using Bo- 
goliubov thcory or its extensions to finite temperature 
0,0. In the sum (g) one should distinguish between 
the collective excitation whose decay is the object of 
the present work and for which we use the notation 
Uosc,v osc ,a OSCl al sc , and the other elementary modes 
which are thermally excited and for which we will use 
the indices i, k already employed in . By keeping in 
(|J) only terms linear in a osc , (aj sc ) and in the produet 
atf.cti (a k ct\), we investigate the desired process where a 
quantum of oscillation fiu is annihilated (created) and 
the z-th excitation is transformed into the k-th. one. The 
same formalism permits to investigate also other decay 
processes where both the excitations i and k are created 
(Beliaev decay of a phonon into two phonons ||) and 
E k + Ei = Huj. In our case we find the following result 
for the rate of energy loss: 



E 



| A ki | 2 8(E k -Ei- hu)(f k - fi) . (5) 



where E = fwjn osc is the energy of the classical oscillation 
{n osc > 1) and 



A ki = 2g J dr %p a [(u* k Vi + v%Vi + u* k Ui)u osc 

+ (v* k v,i + v* k Vi + u* k Ui)v osc \ (6) 

is the relevant matrix element for the process. In deriving 
(|J) we have assumed that at equilibrium the states i, k 
are thermally occupied with the usual Bose faetor f j — 
[exp(j3Ej) — . Introducing the damping rate through 
the relation E = — 2^ E we finally obtain the relevant 
formula 



W7rV|i4 w | 2 S(E k -Ei- hu) 



i k 



9 m) 

dE 



(7) 



where we have further assumed hu> T. Eqs. (|6j[ 
can be further simplificd by caleulating the relevant ma- 
trix elements in semiclassical approximation and can be 
used to caleulate the damping rate of collective excita- 
tions in a trapped Bose gas. This caleulation will be 
the object of a future work. In the following we will use 
them for a systematic diseussion of Landau damping in 
a homogeneous gas where all the ingredients take a sim- 
plified form. In this case i/'o = y/^M) is constant and u, v 
are plane wave funetions: u — U (p) exp(ip • r)/y/V and 



v = V (p) exp(ip • r) /y/V. After some straightforward 
algebra, we can rewrite the rate of the collective exci- 
tation propagating with momentum q in the following 
way: 



7 = 
where 
Bki 



dpi 

(2tt) 3 



\B H \ 2 8{E k -Ei-u) 



dE 



(8) 



2 {[U{E k )V(Ei) + V(E k )V(Ei) 

+ U{E k )U{Ei)]U osc + [V{E k )V(Ei) 
+ V(E k )V(Ei) + U(E k )U(Ei)]V osc } , 



U 2 (E) = l + V 2 (E) = 
U(E)V(E) = -^ 



(E 2 +9 2 nD) 



1/2 



E 



2E 



(9) 



(10) 



and 



E(P) 




gn 



1/2 



2 2 

9 n 



(H) 



is the dispersion law of the elementary excitations in the 
so called Popov approximation JlO|,[llJ , depending on the 
corresponding value of the momentum p. In these equa- 
tions tiq = hq(T) is the condensate density at tempera- 
ture T and the momenta p of the i-th and fc-th excita- 
tions satisfy the condition = p; + q, following from 
momentum conservation. 

A further simplification is obtained if one considers 
the damping for the low frequency excitations, i. e. for 
phonons (q — ► 0). In this case, one has 



V ' \2uj ) 2 \2gnaJ 



1/2 



and 



oscy ' \2u) ) 2 \2gn J 



1/2 



(12) 



(13) 



In the same limit q — > 0, using momentum conserva- 
tion, one can write E k — Ei w v g uj cos(8) / c, where 9 is 
the angle between pi and q, v g is the group velocity of 
the i-th excitation and c is the velocity of sound. After 
integration with respect to 9 and some straightforward 
algebra one finds the useful result 



7 = — CLU 



p 2 dp , c |2 d f (E) 



4irv a 



dE 



where 



C(E) = y/2g[U 2 (E) + V 2 (E) + U(E)V(E) + E 



(14) 



dU 2 {E) 
dE J 
(15) 
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Result (H|) permits to explore explicitly the T- 
dependence of Landau damping. At low T (T <C p, 
where p = gno(0) is the T = value of the chemical 
potential), one finds the result 



7 



21 'w Ujp r , 

"16 J 



(16) 



where p n is the normal part density of the phonon gas 



p A dp 
~6w 2 ' 



df(E) \ 
dE ) 



2w 2 T 4 
45?i 3 c 5 



(17) 



has been first derived by P.Hohenberg and 
01. This eq 



Pn 

Eq. © 

P.Martin |£|]. This eguation can be also obtained fol- 
lowing the method of [13]]. (See ||, Problem to §77.) 

A second important regime occurs at relatively high T 
(T ^> p). Usually at such temperatures the thermody- 
namic quantities are determined by the excitations with 
E ~ T. In the integral ( |l4| ) however the relevant exci- 
tations turn out to have energies E ~ p ~ grto <C T, so 
that one can use the " Rayleigh- Jeans" limit of the dis- 
tribution function, f {E) « T /E. As a consequence the 
resulting dependence of 7 is linear in T and integration 
of (|lj) gives the result 



3w Taq 



(18) 



where we have used the expression c 2 = gn$/m for the 
sound velocity. This regime has been previously investi- 
gated by Szepfalusy and Kondor j^] and more recently by 
Hua Shi [Q . Our result coincides with the one obtained 
in [jTif but the numerical coefhcient slightly differs from 
the one of ||. It is worth noting that result ( p"8| ) docs 
not depend on the value of the condensate fraction tiq. 
The above "high" T regime has been recently employed 
in to provide a quantitative estimate of the width of 
the collective excitations of a dilute Bose gas trapped 
by a harmonic potential. Estimating the momentum q 
of the excitation as Hiu /c where uj is the frequency of 
the excitation and c = (Ti/m)[Awanr,{r = O)] 1 / 2 is the 
value of the sound velocity calculated in the center of 
the trap, the authors of Q| have obtained values for 7 in 
semi-quantitative agreement with the experimental find- 
ings [Q. This agreement suggests that Landau damping 
could represent an important mechanism to explain the 
decay of elementary excitations of trapped Bose gases. 
Notice however that in a non uniform system the elemen- 
tary excitations cannot be described in terms of plane 
wave functions and that Landau damping should be con- 
sequently calculated starting directly from the general 
equation (7), which involves matrbc elements and eigen- 
values relative to the eigenstates of the system. The pres- 
ence of the trapping potential can influence the damping 
mechanism in a deep way. For example the dipole oscil- 
lation does not exhibit any damping in the presence of 
harmonic trapping. 



It is finally interesting to discuss the validity of the 
"high" temperature expansion ([l8|). To this purpose it is 
useful to write the ratio j/ui in the form (see Eq. (14)) 



= (a 3 n (T))^ 2 F(T) 



(19) 



where r = T/mc 2 (T) is a dimensionless variable and the 
function F is given by 



F(r) 



4^ 



2u 2u 2 ' v 



-x/2t\ 



(20) 



with u{x) = (1 + 2; 2 ) 1 



x l 2 . For large r (high temperatures) 
the function F takes the asymptotic value F — > jW Z / 2 t, 
consistent with result ([H]). In Fig. 1 we plot the function 
F versus r. One can see that F approaches the asymp- 
totic linear law rather slowly. For example for tempera- 
tures kT ~ mc 2 (r ~ 1) it differs from it by a faetor 2. 
This suggests that the linear approximation ( |l8| ) should 
be employed with care. In faet in the available traps 0] 
the value of mc 2 , with the sound velocity calculated in 
the center of the trap, corresponds to about 0.4T C [ pd| , 
where T c is the critical temperature for Bose-Einstein 
condensation in the presence of the trapping potential, 
and consequently the condition r 3> 1 is never reached 
below T c . This should be taken into account for quantita- 
tive estimates of the width using the present formalism. 

It is a pleasure to thank S. Giorgini for many useful 
diseussions. 

After completing the present work we have been in- 
formed of a similar analysis of Landau damping JTo| car- 
ried out in dilute Bose gases using semiclassical theory. 
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FIG. 1. Function F plotted as a function of r (full line). 
The linear behavior F — > (3/4)7r 3 / 2 r is also reported (dashed 
line) . 
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